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Time-distance helioseismology is a technique for measuring the time for waves to travel from one point on the solar surface 
to another. These wave travel times are affected by advection by subsurface flows. Inferences of plasma flows based on 
observed travel times depend critically on the ability to accurately model the effects of subsurface flows on time-distance 
measurements. We present a Born- approximation based computation of the sensitivity of time-distance travel times to 
weak, steady, inhomogeneous subsurface flows. Three sensitivity functions are obtained, one for each component of the 
3D vector flow. We show that the depth sensitivity of travel times to horizontally uniform flows is given approximately 
by the kinetic energy density of the oscillation modes which contribute to the travel times. For flows with strong depth 
dependence, the Born approximation can give substantially different results than the ray approximation. 
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1 Introduction 



Time-distance helioseismology (IDuvall et al. 1993b is a 
technique for measuring the time for waves to travel from 
one point on the solar surface to another. Subsurface flows 
advect waves and as a result alter the observed travel times. 
Thus wave travel times can be used as probes of subsurface 
flows (e.g. lKosovichev & DuvaU 1997[[Zhao et al. 2QQ4I) 

One of the key steps in the interpretation of 
travel times is to estimate the effect of subsur- 
face flows on travel times. The ray approximation 
(IKosovichev & DuvaU 19971) , in which travel-time shifts 
are only caused by inhomogeneities located along the 
ray path connecting the observation points, has been 
employed in many time-distance studies of plasma 
flows (e.g. IKosovichev & DuvaU 19"97l IZhao et al. 20011 
IZhao & Kosovichev 2003llZhao et al.^QQ4b . An alternative 
to the ray approximation is the first Born approximation 
(e.g. IBirch & Kosovichev 20001 IGizon & Birch 20021 in 
the context of time-distance helioseimology). The first 
Born approximation takes into account a single scattering 
and thus flows located away from the ray path can affect 
the travel time. 

[Birch & Felder ( 2004)" studied the ranges of validity of 
the Born and ray approximations in a toy problem consisting 
of jets in a homogeneous two-dimensional medium. This 
study showed that there are flow configurations for which 
the Born approximation is valid while the ray approxima- 
tion is not, especially when the transverse size of the jet is 
much smaller than the wavelength. As a result, we would 
like to use the first Born approximation to compute the sen- 



sitivity of travel times to flows for use in time-distance he- 
lioseismology. 

Here we employ the Born approximation approach of 
Gizon & Birch (2002), referred to as GB02 hereafter, to 
compute the sensitivity of travel times to weak, steady, 
three-dimensional subsurface flows in the Sun. We use the 
phenomenological model of Birch et al. (2004), referred to 
as B04 hereafter, to describe how waves are excited and 
damped by convection. We work in Cartesian geometry, 
which is appropriate for waves that travel distances much 
less than the solar radius and also have wavelengths much 
smaller than the solar radius. 

The remainder of this paper is organized as follows. In 
we derive the general expression for the sensitivity of 
travel times to weak flows. In ^we show a few example 
calculations. We compare the Born and ray approximations 
in ^ We conclude in ^ with a summary and a short dis- 
cussion of the implications of the work presented here. 

2 Sensitivity Functions 

In this section we use the Born approximation to obtain the 
linear sensitivity of travel times to weak and steady subsur- 
face flows. We are looking for kernels K = {Kxj Ky^Kz) 
which satisfy 



dT{Xi,X2) 



jjj dr K(r;xi,X2) 



v{r) , (1) 
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where the integration variable r runs over the entire volume 
of the solar model and v = (vx^Vy^Vz) is the flow field. 
Throughout this paper we will denote three-dimensional po- 
sition vectors by r = (x^z) where x = (x^y) is the hor- 
izontal position vector and z is depth. The travel-time dif- 
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ference between surface locations xi and X2 is denoted as 
Sr{x 1^X2) and defined by 



St{xi, X2) = T^{Xi,X2) - T^{x2,Xi) 



(2) 



where r-^{xi^Xj) is the one-way travel time, as defined by 
GB02, from Xi to Xj. 

Following GB02 and B04, we begin by considering 
damped and driven solar oscillations with a displacement 
field ^ that obeys, to lowest order in the flow velocity v. 



(3) 



with the wave equation operator in the absence of flows, Cq 



(e.g. |Lynden-Bell & Qstriker 1967| ), given by 

^oC = PoC - V [7P0 V • C + C • Vpo] (4) 

+(v • Vpo + C • v(Vpo) + Podt (rc) , 

where po, po, and 7 are the background density, pres- 
sure, and ratio of specific heats. We use solar model S 
( [Christensen-Dalsgaard et al. T996). The damping operator 
is r (B04). In equation we have used the Cowling ap- 
proximation and also neglected the variation of the gravita- 
tional acceleration with depth. We have also neglected the 
Coriolis force. The source function S is intended to rep- 
resent the driving of oscillations by near-surface turbulent 
convection. In order to compute the co variance of the wave- 
field ^ we only need the covariance of the source S. We 
choose to use the source covariance described by B04. 

The first-order perturbation to the wave equation intro- 
duced by a flow (e.g. Lynden-Bell & Qstriker 1967) ) is given 
by 

SC^ = 2podtV-V^. (5) 

This term only captures the direct advection of waves by 
the flow V. Any associated changes in the background (non- 
wave) density and sound speed are neglected. The problem 
of the sensitivity of travel-times to changes in sound speed 
has already been addressed by B04. Density perturbations 
could be treated in essentially the same manner. 

The observations used for time-distance helioseismol- 
ogy are typically time series of images of the line-of- sight 
Doppler velocity near the solar surface. For the sake of sim- 
plicity we assume that the line-of-sight is vertical. In ad- 
dition, we approximate the Doppler velocity as the time 
derivative of the displacement at a fixed geometrical height 
Zobs = 200 km. In this case, the observed wavefield, 
(j){x, t), is related to the wave displacement by 



(6) 



where the function denotes the action of the instrument 
point- spread function and any filters applied during the data 
analysis (e.g. phase-speed filters). In the Fourier domain, we 
can write the wavefield as 



(7) 



where k is the horizontal wavevector, uj the angular fre- 
quency, and F{k^uo) represents the filter in the Fourier do- 
main. Throughout this paper we will employ the Fourier 



convention of GB02, for a function / of horizontal position 
and time, we have 



f{x,t) = J dkduj f{k,uj)e'^- 



ik-x—iujt 



(8) 



We will employ the same symbol for functions and their 
Fourier transforms, e.g. f{k^uj) denotes the transform of 
f{x,t). 

In order to compute travel-time kernels we need first to 
obtain Green's functions. We define Green's functions 
as the solutions to 



C^G^ {x, t, z, z') = ej5{x)5{t)5{z - z') , 



(9) 



where with j = z are unit vectors along the coordi- 
nate axes. We use the same boundary conditions as in B04: 
zero Lagrangian pressure at the top of model S and no ver- 
tical motion at a depth of 300 Mm below the photosphere 
(this bottom boundary condition has essentially no effect 
on the computations presented here). The Green's vector 
G^ {x^t^ z^ z') gives the displacement response as a func- 
tion of horizontal position x, time t, and height z, to a delta 
function source in the direction at horizontal position 0, 
time t = 0, and height z\ We use the normal-mode summa- 
tion approximation solution to equation ^ given by B04. 

We now have all of the ingredients to compute travel- 
time kernels using the recipe of equations (26) and (32) 
from GB02. The result for the kernels for velocity is 



Jo 



Kj{r]Xi,X2) =47rRel dujW^^f^{u)Cj{r;x2\xi;u) , 

(10) 

with j = x^y^z. The function IVdiff is the linear sensitivity 
of the travel-time to changes in the cross-covariance and is 
defined in GB02. The sensitivity of the cross-covariance to 
flows in the j direction is denoted by and given by 

Cj{r;x2\xi;uj) = 2{27ryiuj^po{z)m{uj) (11) 
X [IL^ {x — Xi^ z^uj) ■ I{x — X2, z^u) 
+ lL^*{x — X2,z,uj) ■V{x — Xi^z^u)] . 

where r = {x^z). The function m{uj) is the source auto- 
correlation function defined in B04. The vectors I{x) and 
II{x) we define in terms of their horizontal Fourier trans- 
forms: 

V{k,z,uj) = F{-k,uj)Gl{-k,uj,Zohs,z), (12) 
lL^{k,z,Lj) = F{k,uj)H*{k,uj,Zohs)djll*{k,uj,z). (13) 

The Green's vector H is defined as 



H(/c,co',z) = dz'G^{k^uj^ z^ z')\z 



(14) 



where the source depth Zgrc is chosen to be 100 km be- 
low the photosphere (see B04 for a discussion of the source 
model). The general mathematical structure of the kernels 
was explained by GB02. 
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3 Example Calculations 

In this section we show the results of two example calcula- 
tions. The first example is for pi travel-time differences for 
the travel distance A = \\x2 — = 7 Mm obtained us- 
ing a phase- speed filter. The second example is for surface 
gravity wave travel-time differences at a travel distance of 
A = 10 Mm. 

3.1 pi ridge 

For the first example, the filter function F{k^uj) is given as 
the product of three separate filters, 

F{k,uj) = Fi{k,uj)F2{k,uj)0TF{k) . (15) 

The filter Fi removes the /-mode ridge and frequencies be- 
low 1.5 mHz and above 5 mHz. The phase- speed filter, F2, 
is given by 

F2(/c,cj) = e-(^/^-^p)'/'^^' , (16) 

with k = = 12.8 km/s and Svp = 2.6 km/s (this 

is filter 1 from Couvidat et al. 2006). This filter isolates a 
section of the pi ridge. Finally a filter OTF(/c) = e~^^ 
with a = 1.75 Mm is used as a very rough approximation 
for the optical transfer function of the MDl/SOHO high- 
resolution observing mode. 

Figure [T] shows slices through the components of the 
kernel K for travel-time differences, for the pi case. Fig- 
ure [T^ shows a slice through Kx at the photosphere. This 
kernel is symmetric in both x and y. As with the kernels 
shown by GB02 and B04, ellipse- and hyperbola- shaped 
features are visible. The ringing in the horizontal directions 
is a result of the finite band-width of the wavefield. Fig- 
ures [TJ) and [T]: show horizontal slices through the kernels 
Ky and Kz at the photosphere. Different symmetries are 
visible in these slices; Ky is anti- symmetric in both x and y, 
while Kz is symmetric in y and antisymmetric in x. Because 
of these symmetries both Ky and K^ integrate to zero. The 
kernel Kx has a non-zero total integral; spatially uniform 
flows in the x direction cause travel-time differences. 

Figure [T]l shows a vertical slice through Kx dX y = 0. 
Also shown is the ray path corresponding to a frequency of 
4.5 mHz. The maximum height of the ray path is limited 
by the upper turning point (80 km below the photosphere). 
For small travel distances A, such as the example shown 
here, the lower turning point of the ray path is frequency 
dependent. In this case we have chosen to compute a ray at 
4.5 mHz, which is the frequency where the wavefield has 
maximum power (after filtering). The pi mode structure is 
visible in Kx . There is one maximum in sensitivity near the 
photosphere, and another near the lower turning point of the 
mode. This depth dependence will be discussed in later in 
this section. 

A slice through at = is shown in Figure [T^. 
The kernel Kz is largest where the ray is mostly vertical, as 
expected from the ray approximation. 



3.2 Surface gravity waves 

For the second example we compute the sensitivity of /- 
mode travel-time differences to flow, for the case of a travel 
distance A = 10 Mm. In this case the filter function F was 
chosen to be 

Fik.uo) = Fs{k,uj)OTF{k) (17) 

The filter F3 selects only the /-mode and removes all p- 
modes and frequencies below 1.5 mHz and above 5 mHz. 

Figure [21 shows the results of the example /-mode cal- 
culation. The kernels Kx, Ky, and Kz show the same sym- 
metries as in the p-mode case shown in Figure[T] As a result 
the kernel Kx has a non-zero total integral, while Ky and 
Kz both integrate to zero. 

A horizontal slice at the photosphere through Kx is 
shown in Figure [2^. Notice that the Kx kernel is the three- 
dimensional version of the kernel shown by Gizon et al. 
(2000) and in more detail by Jackiewicz et al. (2006). 

Figure [3] shows the depth dependence of the sensitivity 
functions shown in Figures [T] and [21 for the case of hori- 
zontally uniform flows. In both cases, the depth dependence 
is roughly proportional to the kinetic energy density of the 
associated mode. This assumption was suggested for the /- 
mode by Gizon & DuvaU (2000b). 

4 Comparison with Ray Theory 

As described in the introduction, the ray approximation 
JKosovichev & Duvall 1997b has been used to predict the 
travel-time shifts caused by sub- surface flows. In this sec- 
tion we compare the predictions of the ray and Born ap- 
proximations. We will consider here very simple models of 
flows at supergranular scales. We choose to study cylindri- 
cally symmetric flow fields, i;(r, z) = ^^(r, z)z-\-Vr{r^ z)r, 
of the form: 

Vr{r,z) = af{r)hiz) , (18) 
v,{r,z) = ^dr[rfir)]n{z) , (19) 

where r is the distance from the symmetry axis of the flow 
and z is depth. The horizontal variation of the radial flow is 
given by 

/(r) = Jo(fcr)e-'-/^ , (20) 

where L is the decay length of the flow away from the center 
of the cell, and k is the wavenumber associated with the 
radial variation in the flow. The depth variation of the radial 
velocity we choose as 

h{z) = e-(^-^^) Vi^? _ p^-iz-z,r/Di ^ ^21) 

where Zt and are the depths at the top and bottom of the 
cell, Di is the vertical scale of the outflow component of the 
flow, D2 is the vertical scale of the inflow component, and 
the coefficient /3 is chosen so that no vertical flows through 
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Fig. 1 Slices through an example of the sensitivity of a pi travel-time difference to local flows. Panels (a), (b), and (c) 
are horizontal slices, at the photosphere, through the kernels K^, Ky, and Kz respectively. The symmetries of these three 
are different. The kernel Kx is symmetric in both x and y. The kernel Ky is symmetric in x and anti- symmetric for y. The 
kernel is anti- symmetric x and symmetric in y. Because of these symmetries, only Kx has a non-zero total integral, as 
a result the travel-time difference Sr {xi^X2) is not sensitive (at first order) to uniform vertical flows or uniform flows in 
the cross-ray path direction. Panel (d) shows a slice through Kx aty = 0. The heavy black line shows the ray path. Notice 
that the pi mode structure is seen in depth. Panel (e) shows a slice through Kz 3.ty = 0, again with the ray path shown as 
the heavy black line. By symmetry, the kernel Ky is zero at y = 0. In all panels the units are s Mm~^/(km/s). 
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Fig. 2 Slices through the sensitivity of a /-mode travel time to local flows. Panels (a), (b), and (c) are horizontal slices, at 
the photosphere, through the kernels Kx, Ky, and Kz respectively. Panel (d) shows a slice through Kx aty = 0. Panel (e) 
shows a slice through Kz 3.1 y = 0. By symmetry, the kernel Ky is zero aiy = 0. In all panels the units are s Mm~^/(km/s). 
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Fig. 3 Horizontal integrals of the kernels for the pi 
and / cases (solid lines) shown in figures [T] and [21 Also 
shown are scaled kinetic energy densities (dashed lines) for 
the modes at the dominant wavenumber. 

the cell are required by mass conservation. The depth de- 
pendence, n(z), of the vertical flow we then choose so that 
the flow satisfies mass conservation 

dz [pQ{z)n{z)] = po{z)h{z) , (22) 

with the upper boundary condition n{z) = at z = Zt 
where Zt is the top of the cellular flow. 

We consider two models. Both have {L^k) = 
(20 Mm, 0.18 Mm"^). Model A has (zt, ^b, ^i, ^2) = 
(0.2,-8,8,1) Mm and model B has {zt, zi,, Di, D2) = 
(0.2, —2,4,0.5) Mm. In both cases we choose the ampli- 
tude a so that the maximum radial flow speed is 100 m s~^ . 
Model A represents a deep flow that has very little vertical 
variation near the photosphere. Model B represents a shal- 
low flow. The radial and vertical flows for these two models 
are shown in figured 

For both models we compute the travel-time differences 
(5r(iCi, iC2) with (ici, iC2) = (x — A/2, x+A/2)a^, where A 
is the travel distance and x is the unit vector in the x direc- 
tion. For the ray-approximation travel times we use equa- 
tion (15) from Kosovichev & Duvall (1997) with ray paths 
computed according to equation (11) of that paper, for the 
travel distance A = 7 Mm (the same distance as for the ker- 
nels shown in Figure [TJ. We compute Born approximation 
travel-time shifts for the p-mode case described in ^3.11 

Figure [5] shows the ray and Born approximation travel- 
time differences for model A, the deep cellular flow. The 
two approximations give travel-time shifts that are similar 
to within about two seconds. In the Born approximation, 
the contribution of z component of the flow is about 5% of 
the contributions from the x component of the flow and the 
contribution from the y component of the flow contributes 
less then one percent of the signal. We note here that we 
have used, in the ray approximation, only a single ray at 
4.5 mHz (this is the frequency where the wavefield has max- 
imum power). It may be that by taking a weighted aver- 
age of many ray travel times, the ray approximation could 
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Fig. 4 Radial (top panels) and vertical (bottom panels) 
flows for model A (left column) and model B (right col- 
umn). The units of the color scales are m/s. The axis of sym- 
metry of the cellular flow is r = and the photosphere is at 
z = 0. Upflows correspond to positive values of the vertical 
velocity. 
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Fig. 5 Travel-time differences in the ray approximation 
(dashed line) and Born approximation (solid line) for model 
A "deep flow". For this particular flow, the maximum differ- 
ence between the two approximations is 2 s, about 25%. 



be brought into agreement with the Born approximation. A 
study of this procedure is beyond the scope of this simple 
example. 

Figure [6] shows the ray and Born approximation travel- 
time differences for model B, the shallow flow. For this case, 
the ray and Born approximations give substantially differ- 
ent results. The main cause of this difference is that the ray 
approximation is not sensitive to flows that are below the 
lower turning point of the ray, while in the Born approxima- 
tion the sensitivity extends below the lower turning point. In 
this particular example, the strong counter-flow is sensed in 
the Born approximation, but not in the ray approximation. 
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Fig. 6 Comparison between the ray approximation 
(dashed Hne) and Born approximation (soHd Hne) for travel- 
time differences caused by the model B "shallow flow". For 
this flow, the Born approximation and the ray approxima- 
tion give travel-time shifts that differ by up to five seconds. 
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5 Discussion 

We employed the Born approximation to obtain the three- 
dimensional sensitivity of time-distance measurements to 
advection by local flows. In this paper we addressed the im- 
portant question of time-independent mass flows. We have 
shown that for horizontally uniform flows the depth depen- 
dence of the sensitivity of travel times is given approxi- 
mately by the kinetic energy density of the mode which con- 
tributes most to the travel times. 

For simple cylindrically symmetric models of 
supergranulation- scale convection cells we showed 
that the Born and ray approximations can give results 
that are substantially different when the flow varies in the 
depth range just below the lower turning point of the ray. 
This suggests that for inversions of supergranulation- scale 
flows it may be important to use kernels based on the Born 
approximation rather than the ray approximation. 

In the future, a number of improvements could be imple- 
mented: inclusion of modes above the acoustic cutoff fre- 
quency, taking spherical geometry into account, and treat- 
ment of time-dependent flows. 
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